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The long-time behaviors of the time correlation function C v (t) of the shear stress for 
three-dimensional sheared fluids are investigated theoretically. It is found that there are 
the cross-overs in C v (t) from t~ A ^ 2 to t~ 2 for sheared fluids of elastic particles without any 
thermostat, and from t~ 3 ^ 2 to t~ 5 ^ 2 for isothermal sheared fluids including granular fluids. 

§1. Introduction 

The long time behaviors of current autocorrelation functions are important to 
understand the macroscopic properties of fluids P >E}>EKEl i n is known that the ex- 
istence of the long-time tail in the correlation functions leads to the anomalous 
behaviors of the transport coefficients.®'® 

In equilibrium systems, the existence of the long-time tail t~ d / 2 with the time 
t and the spatial dimension d is well recognized. However, the long time behaviors 
of the correlation functions under a steady shear have different feature from those 
at equilibrium.^'® In our previous paper, we find that the velocity autocorrelation 
function C(t) has the cross-over from t~ d l 2 to t~ d in sheared elastic particles without 
thermostat, and C(t) obeys t~( d+2 )/ 2 after the known tail t~ d l 2 in sheared isothermal 
fluids. 81 However, we did not discuss the other correlation functions such as the 
correlation of the shear stress and the related transport coefficients. 

In this paper, thus, we theoretically calculate the correlation function of the shear 
stress. Our theoretical method is based on the classical one developed by Ernst et 
q^EKISJ i n the next section, we will introduce the model we use. In section 3, we will 
present the details of the analysis. In section 4, we will discuss and conclude our 
results. 

§2. Model 

We consider a system consists of N identical smooth and hard spherical particles 
with the mass m and the diameter a in the volume V . The position and the velocity 
of the i-th particle at time t are and Vi(t), respectively. The particles collide 
instantaneously with the restitution constant e which is equal to unity for elastic 
particles or is less than unity for granular particles. When the particle i with the 
velocity Vi collides with the particle j with Vj, the post-collisional velocities v\ and 
v'j are respectively given by v\ = V{ — ^(l + e)(e-Vij)e and v'j = Vj + ^(l + e)(e-Vij)e, 
where e is the unit vector parallel to the relative position of the two colliding particles 
at contact, and V{j = — Vj. 
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Let us assume that the uniform shear flow is stable and its velocity profile is 
given by c a (r) = 'yyS^x, where the Greek suffix a denotes the Cartesian component. 
In this paper, we discuss the following situations: (a) A sheared system of elastic 
particles without any thermostat, (b) a sheared system of elastic particles with the 
velocity rescaling thermostat, and (c) a sheared granular system with the restitution 
constant e < 1. We abbreviate them to the sheared heating (SH), the sheared 
thermostat (ST), and the sheared granular (SG) systems for later discussion. 

We are interested in the correlation function of the shear stress, 

C v (t)= lim l(J v (0)J v (t)), (2-1) 

1/— >oo V 

where Jn(t) is the shear stress at time t. Here, t = is the time when we start 
the measurement. In general, the current Jri{t) consist of the kinetic part and the 
potential part. In this paper, we only consider the contribution from the kinetic part 
of the current. This treatment is correct for dilute gases. For higher density cases, we 
need more sophisticated method to include the contribution from the potential part, 
but the corrections only appear in the prefactor of coefficients at least for elastic 
gases in the equilibrium stateP^ 1 '^ Thus, the currents in Eq. (|2-ip is approximated 
by the kinetic part J^{t) as 

J v (t) ~ tf(t) =Y^mv' ix (t)v' iy (t), (2-2) 

i 

where v[{t) is the peculiar velocity defined as v'^t) = v; L {t) — c(rj(i)). 

§3. Theoretical analysis 

3.1. Hydrodynamic equations 

Following Ernst et a/.,® we approximate the correlation function (|2-ip by 

C v (t)~m 2 n H J dv f o v' 0x v' yfo(v f o ) J dru x (r,t){u y (r,t) - jy), (3-1) 

where fo(v' ) and u(r,t) are the distribution function of the peculiar velocity at the 
initial time t = 0, the velocity field at the position r and the time t, respectively. 
The velocity field u(r,t) obeys the hydrodynamic equations as^ 1 

dtn + V • (nu) = 0, 
d t u + u ■ Vu + (nmY 1 V • II = 0, 
d t T + u-VT + 2(dn)- 1 (n : Vu - k\7 2 T - ^V 2 n) + T( = 0, (3-2) 

where n(r, t) and T(r, t) are the density and the temperature field, respectively. The 
pressure tensor II is given by 

JJij = nTdij - rj (^ViUj + VjUj - VfcUfc^ . (3-3) 
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Note that the bulk viscosity is zero for fluids of dilute hard spheres. £, A, fj,, and 77 
are the cooling rate, the heat conductivity, and the transport coefficient associated 
with the density gradient, and the viscosity, respectively. Here, \x has a finite value 
when e is less than unity but becomes zero in the elastic case. 

Here, the viscosity 77, the heat conductivity A, and the transport coefficient 
associated with the density gradient \x can be non-dimensionalized as 

r/ = 77o7/*, K = KqK* , //=^%*, (3-4) 

n 

where 

f= d{d + 2) r- 

% = aVT, Ko = 2m(d _ 1} aVT (3-5) 

are the viscosity and the heat conductivity in the dilute gasP^ Here, 77*, k*, and n* 
are the constants depend only on e in dilute cases. Here, the explicit form of a is 
given by a = ^^ L r(d/2)n ~ ^/mo~~( d ~ 1 ^ with the gamma function T(x). 

It is obvious that there is the relation £ = for SH. The cooling rate £ becomes 

for ST. On the other hand, the cooling rate is represented by 

C = — C (3-7) 
Vo 

for SG, where £* is the constants proportional to 1 — e 2 . 

3.2. fluctuation of the hydrodynamic fields 

The hydrodynamic equations f|3-2|) have the set of homogeneous solutions as 

n{r,t)=n H , u a (r,t) = jy5 a:X , T(r,t) = T H (t), (3-8) 

where nn = N/V is the homogeneous density. The homogeneous temperature Tn(t) 
is 

T H (t)=T (l + bt)\ (3-9) 

where we introduce the initial temperature To, the characteristic frequency uh = 
uhTh /t]q(Th) , vhq = vh{Tq) and b = 77* 7 2 /(3vho) for SH. On the other hand, Th 
keeps a constant as 

T H (t) = T , (3-10) 
for ST and SG. We can take any value as To for ST, although To depends on 7 as 
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for SG. 

Let us introduce the Fourier transform of any hydrodynamic field f(r,t) is de- 
fined as 

f(k,T) = j d£exp(-fc.£)/(r,r), (3-12) 

where £ is the non-dimensional position defined as £ = Ijj v, with Ijj = 2uh/^ho- 
Let us write the Fourier transport u(k,t) of the peculiar velocity u(r, t) — c(r) 

as 

u(k, t) = fiW (k, t)e^ (fc) + ii( 2 ) (fc, t) e ( 2 ) (fe) + u (3) (fe, t)e^ (fe), (3-13) 

where the unit vectors (fe) depending on k are defined as 

6^ (fe) = k = {k x , ky, k z ) 

e v _ f e (i) . e (!/)) e (l) 

c (2)(fc)- 6 ^ ^ 



|e»-(e« •efe))e( 1 )| 
e ( 3 )(fc) = e « x e< 2 >. (3-14) 

The i-th component of the velocity filed in Eq. (|3- 13[) is expressed as u^(k,t) = 
u(k,t) -eW(fc). 

From the method used in the previous workp^ we can rewrite the velocity field 

as 

fi W(M) = ^(t)j« + {*«<*(*),*) " **>(*(*),*)} 

+ii?(t)^ 1 )(fc(t),o){i?( 1 )(fc(t),t) + ii;( 2 )(fc(t),t)}, 

u( 2 )(fe,i) = fl(t)«W(fc(*),0)£7( 3 )(fe(t),t), 

fi( 3 > (fe, t) = B(t)u® (k(t), 0)£ (4) (k(t), t) + 5(t)^ 2 ) (k(t),0)F(k(t),t), (3-15) 
where we introduce 

fe(*) = (k x , k y - jtk x , kg), (3-16) 
In addition, we also introduce 

^( 2 )(fe,t) = y^e-^^W-T^*)* 2 , 

k(-t) 

E^(k,t) = e-^ k \ (3-18) 
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and 



F(k,t) = M(k(-t))E (3) (k,t) - M(k)E w (k,t) 

with r = rf/2 + + k*, and M(k) = -kk z /(k x k±) taxT 1 (k y / k±) 
Here, we use 



(4) i 



and 



with 



and 



a(t) 



a 2 (t) for SH 

v ai(t) for ST and SG, 

/?(t) = Ai(t) - A 2 (t)k X ky + 4j(t)fc2 



«i(*) = ^-(*y-M-*)*(-*)) 



■— ~ ^gn{^k x ) log 
2jk x 



\k y (-t) - sgn(jk x )k(-t)\ 



\k y - sgn(jk x )\ 



-ai{t), 



+ 



jk x 



Mt) = { 



for SH 
for ST and SG, 



(3-19) 



(3-20) 



(3-21) 



(3-22) 
(3-23) 

(3-24) 



= { 



^ ) iymtf [ 1 + 



2 V * I 7 



9 \» m ) U H° t] [orSr[ 

(i) (^ot) 2 for ST and SG, 



(3-25) 



Mt) = { 



for SH, 



s(^) ( ^°* )3 for ST and SG, 



(3-26) 



where k = yj k x + k% + fc(i) = y/k% + (k y — ■jtk x ) 2 + A; 2 , A; Q = A; a /A;, and 
in Eq. (EFTS!) 

3.3. Result 



Substituting Eq. (|3- 13|) into (|3Tp with the definition of the Fourier transform 
(|3T2p . the correlation function can be written as 



t,i=l,2,3 



(3-27) 
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where we introduce 



C^\t)^m 2 n H l d H J dv' v' 0x v' 0y Mv' ) 



dk 



(2ir) d 

Then, the correlation function can be approximated as 

C v {t)^C^(t) + C^\t) + Cf\t). 



(3-28) 



(3-29) 



('12) ( 21) (31) 

The reason to neglect Cjj (t), (t), and (t) is that these functions decay 
faster than C^ 22 \t) because of the existence of the sound modeP 1 The reason to 
neglect C^ 13 \t) and C^ 23 \t) is that e y 3 \k) is zero. 

In order to obtain C^, (i), C^ 22 \t), and C^ 32 \t), we use Ufc(O) ~ Vo/(noljj), 
and /o(w' ) ~ n(m/27rT ) 3 / 2 exp(-c 2 ){l + a 2 S 2 (c 2 )}, where c = v'^m/2T , S 2 (x) = 
x 2 /2 — 5x/2 + 15/8, and a 2 is the constant depending on the restitution coefficient 
e i4ii2$ Here, we note that a 2 becomes zero when e = 1. Then , C^ 1] (t), C {22 \t), 

(32) 

and C)i (i), are respectively calculated as 



c^\t) = T ^ l + a2) 



lir 3 l d H 



c {22 \t) 



4ir 3 l d H 



C7f)(t)^ r ° 2(1 + a2) 



9 \ 3/2 
9 \ 3/2 
9 \ 3/2 



(3-30) 



where we introduce 



^ u \t) =B 2 {t) j dk 
^ 22 ){t) = B 2 {t) I dk 



k 2 x ky(ky - jtkx) (3( t ) k 2 

k 3 k(-t) 

k 2 ky(k y — jtkx)_ -p(t)k 2 



**(-«) 



^ 2 \t) = B 2 {t) J dUe-J^G(k,t)t 



-/3(t)k 2 



(3-31) 



with 



G(k,t) 



k± 



tan 



1/^2/ 



jtk a 



tan 



(3-32) 



From (EKST]) . V> (11) 0)> ^ 22 \t), and V> (32) (0 obey t" 3/2 for t < Hence, we 



find 



C v (t) oc i 



-3/2 



(3-33) 
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for t < 7 1 . 

However, ^ u \t), ^ 22 \t), and ^ 2 \t) can be evaluated as 

^(ii)( t )~^ c (") 



7' 
86V2 



7 3 t 4 



,(22) 



^ 2 \t) ~ - 



46 l/2 



,(32) 



for t 3> 7 _1 in SH. Here, c^, c^\ and c^ 2 " 1 are the following constants 



(3-34) 



J") 
-if 



.(22) 
'if 



.(22) 
'if 



/ 



dk 



k"xky(ky \^kx) 



(fc2 + fe2)3/2( 2 fc2 + fe 2 + fe 2 _ 2 ^k x k y ) 1 / 2 



^ z k 2 ky(k y \/2kx) c _ 
H - /u^ c 2\/^kxky^ < ^ 



( — ^ — y 



1 + 



-l i ~ V^2^a 



\k z 



tan 



Ik 



tan 



-l / K 



\k z 



(3-35) 



Thus, we find that 



C v (t) oc V (32) (t) oc t 



(3-36) 



for i > 7 _1 in SH. 

On the other hand, ip (11) (t), ip {22) (t), and t/> (22) 0) behave for t > 7- 1 as 



Y>( 22 )(t) 
^ (32) (t)^ 



3\/3 (11 



7 3y;9/2 

3\/3 



7' 



3^9/2 
V3 



(11) 

, 



.(22) 
-7 1 , 



(32) 



Here c^ l \ c^ 2 \ and c£P^ are the constants 



.(22) _ 



'if 



7*5/2 

ts 



(3-37) 



(feg + fc2)3/2( 3A ,2 + fe 2 + fc 2 _ 2^3^^)1/2 



(k "s/^kxky^ 



kg ky ( ky y/3k X ) 



(3A;2 + k 2 + k 2 z - 2^k x k y ) 2 
k ~ 

- ( 

— \~ ky ~\~ k^ ^"sf^kxky 



{k '\/Skxky^ 



{k'^ ^/skxky^ 
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1 + 



tan 



ky \/3kx 

k z 



tan 




Thus, we find that 



C„(t) oc V (32) (t) oct" 5 / 2 



(3-38) 



(3-39) 



for t > 7 _1 in ST and SG. 

Hence, from Eqs. (|3~29|) . (l3^30j) . (13^311) and (13^371) . we find that C v (t) behaves 



as 



t -3/2 £ or t < anci £-2 f or ^ ^ ^, 1 in SH. C v (t) behaves as t~ 3 / 2 for t < 



and r 5 / 2 for t > 7" 1 in ST and SG. 



7" 



§4. Discussion and Conclusion 



In the previous paper,® we find that the exponent of the tail in the velocity 
autocorrelation function C(t) has the crossover from —3/2 to —3 for SH, and —3/2 
to —5/2 for ST and SG. The exponents for ST and SG are common in C v (t) and 
C(t). However, the exponents in C^{t) for SH are different from C(t). This difference 
is originated from the fact that there is B 2 (t) in Eq. (I3-31|) . but the corresponding 
part in C(t) is B{t). 

In conclusion, we have analytically calculated the behaviors of the C n {t) in the 
three dimensional sheared fluids. Based on the method developed by by Ernst et al, 
we find that C„(t) has the long-time tail where the exponent has the cross over from 
-3/2 to -2 for SH and -3/2 to -5/2 for ST and SG. 
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